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Abstract e+e -^tt + tt from the e + e center-of-mass energy W = 

m<p = 1019 MeV down to threshold, i.e. Im-^ < M W7r < 

We present a study of the reaction e+e-^+n'-y at the m ^ ™ s is called radiative return ' The status of the anal y" 
peak with the KLOE detector at the (/.-factory DA$NE. sis and P reliminar y resul * ° n ^ Variant mass spectrum 
This reaction allows us to obtain the cross section for of the two-pion-state are presented. 



1 INTRODUCTION 



the determination of a u 



1.1 Motivation 

The recent surprising results of the g — 2 BNL experi- 
ment [y and the apparent large discrepancy with expec- 
tation, about three times as large as the electroweak cor- 
rection, has renewed the interest in good measurements of 
the cross section e + e~ ^hadrons. The corrections to the 
photon propagator due to low mass hadronic states cannot 
be computed, because QCD is not tractable in this energy 
regime. It is however well known that the correction to the 
muon anomaly, a M , can be related to a(e + e~ — ► hadrons) 
by a dispersion integral. The error on a^(hadr) is in 
fact dominated by the knowledge of the hadronic cross 
section, or equivalently r-decay data. Quoted values for 

([||], using t data) and 
using e + e~ data only). Taking 
theory and experiment differ by 



a M (hadr) are (6924 ± 62)xl0 
(6974±105)xl0- n (jf 
the evaluation of Ref. M 
(426±165)xl0- n . 



1.2 Radiative Return 

We discuss here a method to obtain er(e + e~ — » hadrons), 
which employes the radiative process e + e _ — >hadrons+7, 
where the photon has been radiated by one of the initial 
electrons (positrons). This is referred to as initial state ra- 
diation or ISR [f§], |§]. DA$NE typically operates at a fixed 
center-of-mass energy W = m$. By measuring the radia- 
tive (return) process above, the hadronic cross section be- 
comes accessible over the mass range 2to„. < M n7r < m§. 
The cross section der(e + e~ — ► hadrons)/ds can be ob- 
tained from the measurement of dc(e + e~ — ► hadrons + 
7)/ds, s — W 2 = M^ adr . The relation is given by the 
radiation function H defined by: 

dcr(hadrons + 7) 

s — = er(hadrons) x His, 6 y ) (1) 

ds 

The function H (s, 7 ) (which depends also on the polar 
angle of the photon 7 ) needs to be known to an accuracy 
better than 1%. Radiative corrections have been computed 
so far by different groups up to next-to-leading-order 
for the exclusive hadronic state tt + 7t~ [ph, H, tt/n, 
[ph. We concentrated in the following on the process 
e + e~ — > p7 — > 7r + 7r _ 7, which dominates the cross 
section below 1 GeV. This channel accounts for 62% of the 
hadronic contribution to a^. We also limit our discussion 
to events where the photon is emitted at small angle to 
the beams, 0.934 < cos0 7 < 0.995 (5° < 9 1 < 21°) 
and is thus unobserved. Since the event kinematics is 
only once constrained, large background appears at low 
M nn due to 7r + 7T~7r° and to /i + /i - 7 events, both having 
a cross section larger than 7r + 7r~7. We have also begun 
using large angle photon events, where the kinematics is 
4 times over constrained, three times if we ignore the 75 7 
measurement. The analysis of this channel is in a less 
advanced state but is the best way for studying the low 
dipion mass region, which is however very important for 



The radiative return method has one important ad- 
vantage over an energy scan. The systematics of the 
measurement (e.g. normalization, beam energy) have to 
be taken into account only once while for the energy scan 
they have to be known for each energy point. 

1.3 Final state radiation 

Final state radiation from pions is an intrinsic background, 
different from background processes such as muon pairs, 
machine background etc. and must be distinguished from 
ISR. It is in general possible to correctly extract the rele- 
vant 7r + 7r~7 contribution from the FSR background with a 
fit to the data of the cross section for both processes, after 
proper folding of the detector response. In the following 
we suppress FSR events by kinematic cuts, which reduce 
this background to an acceptable limit. 
Monte Carlo simulations [nQ| show, that cutting on E 1 = 
(m0 — E + — E_) and 7 , the energy and polar angle of 
the photon in the e + e~ system, effectively suppresses FSR, 
while most of the ISR events are retained. ISR events are 
peaked at large cos 7 , while FSR photons tend to be emit- 
ted along the pion direction ex sin 2 ^. 

A cut on E 1 additionally suppresses FSR due to the fact 
that the decay via the p resonance (i.e. ISR) leads to an 
enhancement of the photon energy spectrum at s=s 220 MeV 
at which value FSR is very suppressed. The cuts used in the 
analysis given below reduce FSR to less than 1%: 



0.934 < cos 7 < 0.996 
E 7 > 10 MeV 
55° < 9* < 125° 
p^ > 200 MeV/c 



(2) 
(3) 
(4) 
(5) 



The effective cross section with these acceptance cuts is 
4.2 nb. Note that, from the experimental point of view, 
it would be much easier to perform an analysis with no 
lower boundary for 6 1 (0° < 7 < 21°). This would 
increase the cross section by a factor ~ 4 and would 
eliminate a detector resolution effect due to the lower 
angular cut. However, this is not correct in the current 
version of the generator, which diverges for 6> 7 — ► 0°. A 
new calculation with the NLO QED corrections taken into 
account has recently been published[0| and a new version 
of the generator will soon be available. This will allow 
to remove the lower angular cut improving the statistical 
error and removing a possible source of systematics. 

The description of FSR is model dependent and in 
the actual version of Ref. [g] a point like pion is assumed. 
The model dependence can be tested for 7r + 7r _ 7 events by 
looking at the charge asymmetry of the produced pions: 



A(6„) 



N* + (8*) + N*- (0^) 



(6) 
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Figure 1: The charge asymmetry (formula (|6|)) is shown 
for 2 angular regions of the photon polar angle: (a) 60° < 
6^ < 120°, where we see a sizable effect of the charge 
asymmetry due to larger FSR and (b) 5° < 7 < 21°, 
where FSR and hence the charge asymmetry are small. 

This charge asymmetry arises from the interference be- 
tween ISR and FSR and is therefore linear in the FSR am- 
plitude. We measured the charge asymmetry for (a) large 
photon angles (60° < 7 < 120°), where FSR is large, and 
(b) for small photon angles (5° < 7 < 21°). The results 
are illustrated in figure (Q) and show very good agreement 
between data and Monte Carlo, indicating that the point 
like model describes well the process of FSR within the 
error bars. 

1.4 Studying the Pion Form Factor at threshold 

For the selection cuts described above (small angle region), 
the ISR cross section around the p peak is highest, but the 
M%^ region at threshold is suppressed by the kinematical 
cuts. To avoid this, a second angular region has been de- 
fined in which we can explore the invariant mass spectrum 
down to s/Wj^ ~ 300 MeV: 



As discussed above, FSR is in average much higher at 
large values of 7 , but does not play a role anymore for 
small values of M% w , which we actually want to study with 
these new cuts. 



2 EVENT SELECTION 

In this chapter we present the event selection for the mea- 
surement of the 7r + 7r~7 final state. After a very brief de- 
scription of the KLOE detector, the selection algorithm for 
this signal is presented. 



2. 1 The KLOE detector 



< cos0 7 < 0.5 
E 1 > 10 MeV 



(7) 
(8) 



KLOE [ |11[ ] is a typical e + e~ multiple purpose detector 
with cylindrical geometry, consisting of a large helium 
based drift chamber (DC, Jl2j]), surrounded by an electro- 
magnetic calorimeter (EmC, [|l3|]) and a superconducting 
magnet (B = 0.6 T). The detector has been designed 
for the measurement of CP violation in the neutral kaon 
system, i.e. for precise detection of the decay products 
of Ks and K^. These are low momenta charged tracks 
(7r ± , |U, , e* with a momentum range from 150 MeV/c to 
270 MeV/c) and low energy photons (down to 20 MeV). 
The DC dimensions (3.3 m length, 2 m radius), the drift 
cell shapes (2x2 cm 2 cells for the inner 12 layers, 3x3 
cm 2 cells for the outer 46 layers) and the choice of the 
gas mixture (90% Helium, 10% Isobutane; X = 900 
m) had to be optimized for the requirements prevailing 
at a factory. The KLOE design results in a very good 
momentum resolution: a p± /p± < 0.3% at high tracking 
efficiencies (> 99%). 

The EmC is made of a matrix of scintillating fibres 
embedded in lead, which guarantees a good energy re- 
solution cteJE = 5.7%/y/E(GeV) and excellent timing 
resolution a t = 57ps/y/E(GeV) 50 ps. The EmC 
consists of a barrel and two endcaps which are surrounding 
the cylindrical DC; this gives a hermetic coverage of the 
solid angle (98%). However, the acceptance of the EmC 
below sa 20° is reduced due to the presence of quadrupole 
magnets close to the interaction point and does not allow 
to measure e.g. the photon of tt + it~j events with low 7 
angles (as required for FSR suppression). 

It will be shown in the following, that an efficient se- 
lection of the 7r + 7r~7 signal is possible, without requiring 
an explicit photon detection. The relatively simple sig- 
nature of the signal (2 high momentum tracks from the 
interaction point) and the good momentum resolution of 
the KLOE tracking detector allow us to perform such a 
selection. 
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Figure 2: A likelihood method has been developed to separate electrons from pions. The kinematic variable track mass 
before and after the application of this method is shown. 7r + 7r~7 events are peaked at MTrack — rn^, radiative Bhabha 
events at much smaller values. Other background channels are also visible (/i + /i~7 and ir + Tr~ir°). 



2.2 Selection Algorithm 

The 7r + 7r~7 events are selected using the following 4 
steps. The selection is based on the measurement of the 
charged pion tracks by the DC, while the photon is not 
required to be detected in the EmC. Calorimeter informa- 
tion is however used for the 7r/e-separation (likelihood 
method). 



i) Charged vertex in DC. We require 1 vertex in the 
DC with 2 associated charged tracks close to the in- 
teraction point: y/{x"y + y v ) < 8 cm and \zy\ < 15 
cm. 

ii) Likelihood Method for ir/e - separation. A fraction 
of radiative Bhabha events e + e~7 enters the kinematical 
selection (see next selection cut), giving rise to a non neg- 
ligible background. In order to reject those events, a like- 
lihood method has been worked out for an effective ir/e- 
separation. The method is based on the shape and energy 
deposition of the EmC clusters produced by the charged 
tracks and has been developed using independent control 
samples for the pion information (7r + 7r~7r° events) and 
for the electron information (e + e~j events). 98% of all 
7r + 7r~7 events are selected if at least one of the two tracks 
has been identified as a pion by the likelihood method. In 
figure (g) the effect of the likelihood method is demon- 
strated in the track mass (MTrack) distribution!]. 7r+7r_ 7 
events are peaked at MTrack — m-K, radiative Bhabha 
events at smaller values. 

iii) Kinematic Cut: 130.2 MeV < M Track < 
149.0 MeV. We perform a ±9.6 MeV (±2ct) cut on 
the kinematical variable Mxracki which is peaked at 



The Mx ra ck variable is obtained applying the 4-momentum- 
conservation and resolving it for the particle mass. 



"7 events (see figure (0)). Back- 
ground of 7r + 7T~7r° events (with MTrack mostly > 150 
MeV), ^i + /i _ 7 (with MTrack peaked at m,u) and the bulk 
part of radiative Bhabha events e + e~7 (Mxrack < 100 
MeV) are mostly rejected (see subchapter 3.1 on back- 
ground). This cut assumes a ir + ir~j final state. 

iv) Acceptance Cuts. The missing momentum of the 2 
accepted charged tracks is calculated and associated with 
the photon under the assumption of a 7r + 7r~7 event: p 1 = 
P0 — p+ — pi n. The acceptance cuts of formulae (0)-(|l) 
are then applied. 



3 EVENT ANALYSIS 

The 7r + 7r~7 cross section measurement contains the fol- 
lowing terms: 



da h 



adrons+7 



dM^ 



dNp hs - dJV; 

dM 2 



Bk s 



1 



esci£Acc / Cdt 



(11) 



It requires the study of the various background channels 
(iVekg) see following subchapter), the selection efficiencies 
(esei) and the systematic effects due to the acceptance cuts 
applied (eacc)- Finally the counting rate measurement has 
to be normalized to the integrated luminosity J Cdt in or- 
der to achieve a cross section measurement. All these terms 
will be obtained from data for the final analysis. 
Preliminary results concerning esd are presented in the fol- 
lowing subchapter. The detector resolution and the system- 
atic effects, arising from there have been studied in detail 
with MC [[[(J . No limititations have been found for a high 
precision measurement at the percent level. 



pi is the c 
crossing angle. 



momentum in the laboratory system due to the beam 



3. 1 Background 

Radiative Bhabha events are mostly suppressed by the 
Mxrack cut (chapter 2.2 (ii)). The remaining background 
due to events with a high value for Mxrack an d due to 
electrons, which have not been rejected by the likelihood 
method, is peaked at large M%^ (above 0.7 GeV 2 ) and 
corresponds to a contamination below the percent level in 
this M 2 „. region. 

/i + /i~7 events are not efficiently rejected by the likelihood 
method, because they release a pattern in the EmC with a 
signature similar to pions. After the cuts of formulae (ph - 
(Bb and after the Mxrack cut , the remaining n + [i~ cross 
section is low (as 1CP 2 nb), such that we expect only a 
small contamination (< 1% in the high M % v region > 0.7 
GeV 2 ). 



An important background for our signal is the de- 
cay (j> -> Tr+ir-Tr (B.R. 15.5%, <+ 7r - 7r o ~ 500 
nb). 7T + 7r~7r° and ir + ir~j events are separated in the 
KLOE standard reconstruction scheme by a cut in the 
2-dimensional plane Mxrack — M% v . At small M% n , the 
Mxrack values for the 2 channels are very similar and a 
part of the Tc + Tr~ir° events appear as background. The 
Mxrack cut an d the acceptance cut of formula (^) [] reject 
a big part of these events. 

In order to estimate from data the remaining contamina- 
tion, we modified the standard cut in the Mrrack ~ M^ 
plane by expanding the tt + tt^tt selection region. We 
see then the tail of ir + ir~ir° events entering the Mrrack 
selection interval. We perform this study in bins of M 2 ^. 
The 7r + 7T~7r° background is negligible in most of the M 2 OT 
region and gives only a contamination at the lower end 
of the spectrum between 0.3 GeV 2 and 0.4 GeV 2 . The 
effective 7r + 7r~7r° cross section after all the selection steps 
is < 0.01 nb. It increases considerably if we select 7r + 7r~7 
events at larger polar angles of the photonn 

3.2 Selection Efficiencies 

We shortly summarize the various efficiencies which con- 
tribute to the total selection efficiency: 

• Trigger: The trigger efficiency has 2 contributions: 
the probability of a 7r + 7r~7 event to be recognized 
by the KLOE trigger (between 95% to 99% depend- 
ing on M 2 „.) and the inefficiency which arises from 
a trigger hardware veto for the filtering of cosmic ray 
events. The second contribution causes an inefficiency 
for 7r+7r _ 7 events at large M 2 T (30% at 1 GeV 2 ) 
which decreases with lower M% v (fully efficient at 
« 0.7 GeV 2 ). These values have been obtained from 



3 The pion tracks have on the average a lower momentum as ■n + -n~'y 
events. 

4 This behaviour can be easily explained by the missing momentum of 
7T+7T - 7r° events (in this case the ir°), which is peaked at large angles. 



data by looking at the individual probabilities for 7r + 
and 7r~ to fire 1 trigger sector and 1 cosmic veto sec- 
tor. 

• Reconstruction Filter: A software filter is imple- 
mented in the KLOE reconstruction program for the 
filtering of non-collider physics events, like e.g. ma- 
chine background and cosmic ray events. The ineffi- 
ciency for 7r + 7r~7 events caused by this filter is m 2% 
(taken from MC). 

• Event Selection (see subchapter 2.2): The DC vertex 
efficiency (s=s 95%) is obtained from Bhabha events, 
selected without requiring DC information. The ef- 
ficiency due to the likelihood selection is « 98% 
and is evaluated from data during the construction 
of the likelihood method. The efficiency due to the 
Mxrack cut (~ 90%) is evaluated from Monte Carlo 
at present. 

3.3 Luminosity Measurement 

The DA$NE luminosity is measured by KLOE using large 
angle Bhabha (LAB) events. The effective Bhabha cross 
section at large angles (55° < #+,- < 125°) is still as 
high as 425 nb. The number of LAB candidates Nlab 
are counted and normalized to the effective Bhabha cross 
section, obtained from Monte Carlo: 

^= ^#^ -(l-W (12) 

<J LABV j +-) 

Hence, the precision of this measurement depends on: 
(i) the theoretical knowledge of the Bhabha scattering 
process including radiative corrections; (ii) the simulation 
of the process by the detector simulation program. 

For the theory part we are using 2 independent Bhabha 
event generators: the Berends-Kleiss [Q generator, 
modified for DA$NE in [JL5J] and the BABAYAGA 
generator[|lq]. 

We use a selection algorithm for LAB events with a 
reduced number of cuts, for which we expect a very good 
description by the KLOE detector simulation program. 
The acceptance region for the electron and positron polar 
angle (55° < + ._ < 125°) is measured by the EmC 
clusters produced by these tracks, while the energy mea- 
surement (25+ > 400 MeV) is performed by the high 
resolution drift chamber. Taking the actual detector reso- 
lutions, we expect the systematic errors arising from these 
cuts to be well below 1%. Background from /i + /u — (7), 
7r + 71-^(7) and 7T + 7r~7r° is below 1% and can be easily 
subtracted. All the selection efficiencies concerning the 
LAB measurement (Trigger, EmC cluster, DC tracking) 
are above 98% and are well reproduced by the detector 
simulation. 

As a goal we expect to measure the DA$NE luminosity 
at the level of 1%. The very good agreement of the 
experimental distributions (0+,-, £+,-) with the existing 



event generators and a cross check with an independent 
luminosity counter based on e + e~ — ► 77(7), indicate a 
good accuracy. However, more systematic checks (e.g. 
the effect of a varying beam energy and of a non centered 
beam interaction point) are still to be done. 
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Figure 3: The differential ir + ir~j cross section as a func- 
tion of the pion invariant mass. The solid line is the predic- 
tion of the theoretical event generator EVA of Ref . [ra] . 



3.4 Comparison with Monte Carlo 

We analyzed a data sample of 16.4 pb _1 and present in fig- 
ure [3] the preliminary result for the dipion invariant mass 
spectrum M 2 ^ (small angle region) and in figure |4J the pion 
form factor IF^] 2 . Therefore we applied formula (p and 
write for awrons = na 2 {3l/{3M 2 „) x \F V (M^)\ 2 . In 
plot |] both the small angle (formulae (gj)-(|5|)) and the large 
angle events (formulae (0)-(|lO|)) are plotted. The solid line 
is a Gounaris-Sakurai parametrization with a fixed set of 
parameters and allows a qualitative comparison with data. 
No attempt has been made to subtract the residual contri- 
bution of <fi — ► ir + ir~ir° events which are present at low 
M 2 

7T7T 

The statistical error of the data points in the p peak re- 
gion is rs 2%. The actual version of the event generator 
has a systematic uncertainty of the same size. By com- 
paring the data and MC distributions we conclude that the 
accuracy of our cross section measurement is at the level of 
a few percent, which will be considerably improved with 
the ongoing analysis. 



the radiative return method. The data sampleR, correspond- 
ing to an integrated luminosity of 16.4 pb _1 shows a good 
agreement with the "standard" parametrization for the p 
peak. It also shows the possibility for KLOE of explor- 
ing the hadronic cross section at threshold. We conclude 
that the experimental understanding of efficiencies, back- 
ground and luminosity are well under control. We further 
conclude that the radiative return is a competitive method 
to measure hadronic cross sections. We note again the ad- 
vantage of the radiative return. Systematic errors, like lu- 
minosity and beam energy, enter only globally and do not 
have to be known for individual energy points. 

For the future we plan to refine the analysis and to in- 
clude very small angle photons. A comparison with MC 
will be possible in this case with the new next-to-leading- 
order event generator [J7]] and will improve the precision of 
the measurement due to a better systematic control of the 
fiducial volume. Moreover we are studying the possibility 
to enlarge the acceptance region for 0+,-, which increases 
the kinematical acceptance of events at low M 2 ^ , <0.3 
GeV 2 . In order to improve on the accuracy of a M (hadr) 
and to be competitive with results coming from the CMD- 
2 experiment in Novosibirsk [|l7[, a final precision for this 
measurement below 1% is needed. A statistical error on 
this level can be obtained with an integrated luminosity 
of ^200 pb _1 which could be collected before year end. 
DA$NE has recently achieved a luminosity £>1.5 pb _1 

We are also studying how to measure d<r(e + e~ — > 
7r + 7r - + nj))/dM 2 n , which simply means chosing events 
which satisfy appropriate kinematical constraints to reject 
unwanted backgrounds and then apply precise tranverse 
momentum and energy unbalance for proper inclusion of 
soft photons, to facilitate computation of the physical cross 
section. The radiative corrections for such a measurement 
are obtainable with high precision (< 0.5%)|j. 

At low energy, Af W7r < 600 MeV, the p + p~"[ 
cross section becomes larger than the dipion cross sec- 
tion. It therefore becomes possibile to measure the ratio 
er Tr + 7r - 7 /cr Al + M - 7 in which initial state radiation and vac- 
uum polarization effects cancel. Final state radiation cor- 
rections are still necessary. On the experimental side one 
needs only to estimate the difference in the efficiency for 
the two channels rather then the entire efficiency. 
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4 SUMMARY AND OUTLOOK 
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We presented in this paper a preliminary measurement 
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< 1 GeV 2 , using 



5 This corresponds to about one half of the full data set which KLOE 
has taken in 2000. 
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Figure 4: The pion form factor as a function of the y/s of the hadronic system for small and large angles. The line is a 
Gounaris-Sakurai parametrization with fixed parameters. No fit has been attempted. 
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